#Metview Macro

t_fc24 = read("t_fc24.grib")

t_fc96 = read("t_fc96.grib")

diff = t_fc96 - t_fc24

pos = mcont(
    legend                         : "on",
    contour_level_selection_type   : "level_list",
    contour_max_level              : 10,
    contour_min_level              : 0.5,
    contour_level_list             : [0.5,1,2,4,10],
    contour_shade                  : "on",
    contour_shade_method           : "area_fill",
    contour_shade_max_level_colour : "red",
    contour_shade_min_level_colour : "orange_yellow",
    contour_shade_colour_direction : "clockwise"
    )

neg = mcont(
    legend                         : "on",
    contour_level_selection_type   : "level_list",
    contour_max_level              : -0.5,
    contour_min_level              : -10,
    contour_shade_max_level        : -0.5,
    contour_shade_min_level        : -10,
    contour_level_list             : [-10,-4,-2,-1,-0.5],
    contour_shade                  : "on",
    contour_shade_method           : "area_fill",
    contour_shade_max_level_colour : "greenish_yellow",
    contour_shade_min_level_colour : "blue",
    contour_shade_colour_direction : "clockwise"
    )

xs_europe = mxsectview(
    line : [55,-6,43,16]
    )

plot(xs_europe,diff,neg,pos)